Bloch-like oscillations in a one-dimensional lattice with long-range correlated disorder 
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We study the dynamics of an electron subjected to a uniform electric field within a tight-binding 
model with long-range-correlated diagonal disorder. The random distribution of site energies is 
assumed to have a power spectrum S(k) ~ l/k a with a > 0. Moura and Lyra [Phys. Rev. Lett. 81, 
3735 (1998)] predicted that this model supports a phase of delocalized states at the band center, 
separated from localized states by two mobility edges, provided a > 2. We find clear signatures 
of Bloch-like oscillations of an initial Gaussian wave packet between the two mobility edges and 
determine the bandwidth of extended states, in perfect agreement with the zero-field prediction. 

PACS numbers: 78.30.Ly; 71.30.+h; 73.20.Jc; 72.15.Rn 



The single-parameter scaling hypothesis predicts local- 
ization of a single quasiparticle in one (ID) and two di- 
mensions with time-reversal symmetry, independently of 
the disorder strength present in the system [lj . There ex- 
ist, however, low-dimensional systems that do not obey 
the single-parameter scaling framework. Thus, the ab- 
sence of Anderson localization in the presence of spatial 
short-range correlations in disorder 0,13 was put forward 
to explain transport properties of semiconductor super- 
lattices with intentional correlated disorder £|. Further, 
it was demonstrated that long-range correlated diago- 
nal 0,01 and off-diagonal Q disorder also acts towards 
delocalization of ID quasiparticle states. Furthermore, 
long-range correlations can result in the emergence of 
a phase of extended states in the thermodynamic limit. 
This phase appears at the band center and is separated 
from localized states by two mobility edges 5] . This the- 
oretical prediction was experimentally validated by mea- 
suring microwave transmission spectra of a single-mode 
waveguide with inserted correlated scatterers [8j. In the 
case of short-range correlated disorder the delocalized 
phase does not appear: the number of delocalized states 
increases proportionally to the square root of the system 
size, and thus this phase has zero measure in the ther- 
modynamic limit. 

In this Letter we focus on the dynamical properties of 
an electron in a system with long-range correlated diago- 
nal disorder. Interplay between the delocalization effect, 
preserved by the long-range correlated disorder, and the 
dynamic localization, caused by an electric field acting 
on the system, is of our interest. We compute the behav- 
ior of an initial Gaussian wave packet in the presence of 
a uniform electric field solving numerically the ID time- 
dependent Schrodinger equation for the complete Hamil- 
tonian. We found clear signatures of Bloch-like oscilla- 
tions of the wave packet between the two mobility 



edges of the delocalized phase of states. The amplitude 
of the oscillatory motion of the centroid allows us to de- 
termine the bandwidth of the delocalized phase. 

We consider a tight-binding Hamiltonian with long- 
range-correlated diagonal disorder and an external dc 
electric field on a regular ID open lattice of spacing a 
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where |n) is a Wannier state localized at site n with en- 
ergy e n , J- is the external uniform electric field, and N 
is even. The intersite coupling is restricted to nearest- 
neighbors and assumed to be uniform over the entire lat- 
tice with J > 0. In terms of the Wannier amplitudes 
ip n (t) = (n\^>(t)), the Schrodinger equation reads 



iip n = ( e „ - Fn)ip n - ip n 
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where we introduced the dimcnsionless parameters e n = 
s n /J, F = eTa/J, and time is expressed in units of h/J. 

The source of disorder is the stochastic fluctuations 
of energies e n , which we are going to consider as being 
long-range correlated. One of the simplest ways to nu- 
merically generate a power-law correlated sequence of on- 
site potentials e„ is to write its Fourier decomposition as 
follows @ 
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Here, (j) k are N/2 independent random phases uniformly 
distributed within the interval [0, 2tt], and C{a) is a nor- 
malization constant. We will normalize the energy se- 
quence to have (e n ) = and (e^) = 1, where (...) in- 
dicates average over the random phases The long- 
range nature of the potential correlations results from the 
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power-law dependence of the amplitudes on the wave- 
vector characterizing each Fourier component. Several 
stochastic processes in nature are known to generate 
long-range correlated random sequences which have no 
characteristic scale, for example, in the nucleotide se- 
quence of DNA molecules . The relevance of the un- 
derlying long-range correlations for the electronic trans- 
port in DNA has been recently discussed in Ref. . Fur- 
thermore, interface roughness appearing during growth 
often displays height-height correlations with power-law 
spectra [l4j ; thus, the subsequent random potential aris- 
ing from the rough interface would be long-range cor- 
related. Recently, transport properties of systems with 
long-range correlated disorder was explored, both theo- 
retically and experimentally, in the desi gn o f devices for 
filtering of electrical and optical signals [l5| . 

Having introduced the model of disorder, we can solve 
numerically J5J by means of an implicit integration al- 
gorithm [13 ■ The initial condition is set to a Gaussian 
wave packet of width a and centered at no = N/2 

i/> n (0) = A(a) exp [-(n - n )2/4a2] , (4) 

A(cr) being the normalization factor. Once Eq. J2J is 
solved for the initial condition Q , we compute the mean 
position of the wave packet (centroid) 
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It is to be noticed that all eigenstates of the Hamilto- 
nian contribute to x(t). 

It is well known that a uniform electric field applied 
to a periodic lattice causes the dynamic localization of 
the initially extended states of an electron and gives rise 
to an oscillatory behavior of the electron wave packet, 
the so-called Bloch oscillations ||. The length of the 
segment over which the electron oscillates (twice the 
amplitude) and the period of these oscillations are es- 
timated from semiclassical arguments to be (in dimcn- 
sionless units) Lp = W/F and tb = 2tt/F respectively 
(e.g., see Ref. ^^), where W is the width of the Bloch 
band in units of the coupling integral J (W = 4 in our 
case). Notice that this approach requires slowly varying 
wave functions to be valid, implying that Lp > 1. Simi- 
lar condition was pointed out in Ref. The question 
we aim to clarify is to what extent will the correspond- 
ing phenomenology be valid if long-range correlated dis- 
order, that allows for a delocalized phase, is present. In 
other words, whether the biased phase of extended states 
bounded by two mobility edges behaves similarly to a bi- 
ased Bloch band. If so, it provides us with a method to 
measure the energy width of the delocalized phase W 
from the relationship L' F — W' /F, where L' F is twice 
the amplitude of the oscillatory motion of the centroid. 
Below, we present a numerical proof of this conjecture. 

In Fig. n we show the asymptotic dynamics [t ^> Tb) 
of the Wannier amplitudes tp n {t) of a initial Gaussian 



wave packet ^ n {t) (a = 10 at t = and F — 0.1) cal- 
culated for two values of the exponent a, 0.5 (the phase 
of delocalized states is absent in the unbiased system) 
and 3.0 (the phase of delocalized states is present in the 
unbiased system). At a = 0.5 [see Fig.^i)], the disorder 
is almost uncorrelated and the system is similar to the 
standard Anderson model, with no signatures of Bloch 
oscillations. On the contrary, the plot in Fig. QJ)) cal- 
culated for a = 3.0 shows an oscillating in time pattern. 
This result suggests that Bloch oscillations can take place 
even in the presence of disorder. 




FIG. 1: Asymptotic dynamics (t 2> tb) of the Wannier am- 
plitudes ip n (t) of a biased wave packet (a = 10 at t = and 
F = 0.1) in a lattice with TV = 300 sites, for a) a = 0.5 and 
b) a — 3.0. Bright and dark regions indicate higher and lower 
probability amplitudes, respectively. 

The time-domain evolution of the centroid, x(t), pro- 
vides more detailed information about what is happen- 
ing. In Fig. |21 one can see that there are no signatures 
of Bloch oscillations for a — 0.5 (dotted line), at least 
for a moderate field amplitude. Oscillations, which are 
present at the very beginning, achieve in a short time a 
weakly fluctuating (stationary in average) value. This is 
a clear indication of the existence of disorder-induced de- 
coherence effects, as is expected for the localized regime. 
On the contrary, at a = 3.0 the centroid displays an 
oscillatory, amplitude-modulated pattern after an initial 
transient, where the amplitude of the oscillations is re- 
duced (solid line). This oscillatory displacement of the 
wave packet is not accompanied by in-phase oscillation 
of its width. In other words, breather modes are absent 
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for the specified initial condition. 
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FIG. 2: Time-domain dynamics of the centroid of a biased 
wave packet (a — 10 at t = and F — 0.25) in a lattice with 
N = 500 sites. Dotted lines correspond to a — 0.5 and solid 
lines to a — 3.0. 



Valuable information can be extracted from a detailed 
inspection of the centroid dynamics. First, we notice that 
the period of the oscillations is a well defined quantity, 
as is seen from Fig. [3J Furthermore, within the numer- 
ical uncertainty, its value is equal to tb = 2tt/F, being 
the period of Bloch oscillations in a homogeneous lat- 
tice. At the initial stage of the motion, the spatial region 
within which the wave packet oscillates is roughly given 
by 4/F, as for an ideal Bloch band. Disorder results 
in a relatively fast decrease of the oscillation amplitude. 
After this transient stage, oscillations are not damped 
further and remains amplitude-modulated with an enve- 
lope which depends on the realization of disorder. Dotted 
lines in Fig. [3] bound the spatial region within which the 
wave packet oscillates for a long time. On average, the 
width of this region L' F is found to be L' F ~ W /F, where 
W' is independent of the applied field F. From the data 
in Fig. |31 we obtain W' ~ 1. This value agrees remark- 
ably well with the width of the band of extended states 
reported in Ref. ,5] for a = 3.0. Thus, we arrive at the 
main conclusion of this work, namely there exist clear 
signatures of Bloch-like oscillations of a biased Gaussian 
wave packet between the two mobility edges. 

To provide further confirmation of this claim, we calcu- 
lated numerically the Fourier transform of the centroid, 
x(u>), as shown in Fig.HJ Results were obtained by aver- 
aging over 1000 realization of disorder. For a — 0.5, the 
Fourier transform x(w) is rather broad, suggesting that 
x(t) is similar to a white noise signal. On the contrary, for 
a = 3.0 the Fourier transform x(u>) shows a well-defined, 
narrow peak at about u> = F, despite averaging. The 
variance of x(u>) , cr^ > displays a rather different trend be- 
low and above a = 2.0, as seen in the inset of Fig.^J For 
a < 2.0 the Fourier transform is broad, but for a > 2.0 
the width of the Fourier transform is rather small and 
independent of the electric field. The onset for the ap- 
pearance of Bloch oscillations a = 2 (see the inset of 
Fig. 0J) is in excellent agreement with the value obtained 
in Ref. 5] for the occurrence of the phase of extended 
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FIG. 3: Time-domain dynamics of the centroid of a biased 
wave packet (a = 10 at t = 0) in a lattice with iV = 1000 
when a = 3.0, for a) F = 0.125 and b) F = 0.5. Dotted 
lines bound a region of size L' F within which the wave packet 
oscillates for a long time. 




FIG. 4: Fourier transform of the centroid of a biased wave 
packet (a = 10 at t = and F — 0.5) in a lattice with 
N — 500 sites, for a = 0.5 (upper panel and notice the mag- 
nification factor) and a — 3.0 (lower panel). The inset shows 
the variance ct5 of x(ui) as a function of a for different val- 
ues of the applied electric field. Results were obtained by 
averaging over 1000 realizations of the disorder. 



states. 

Before concluding, some words concerning the influ- 
ence of the width of the initial wave packet are in order. 
The above results for F ^ were obtained with a rel- 
atively broad wave packet (a 3> 1) with initial velocity 
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FIG. 5: Time-domain dynamics of the centroid of a biased 
wave packet, with initial condition ipn(0) — 8 nnQ1 in a lattice 
with N = 1000 sites when a = 3.0, for F = 0.125 



equal to zero. Consequently, its Fourier components span 
a very narrow set of states at the bottom of the band. 
We performed also a study for initial moving wave pack- 
ets located initially deep in the band. In short, in the 
very beginning the wave packet explores the entire band 
in either case and oscillates within a spatial region of size 
4/F. On increasing time, the amplitude of x(t) is always 
reduced to a region of size L' F — W'/F, where W is 
the width of the delocalized phase. This picture changes 
slightly for narrower initial wave packets, for which the 
large-amplitude transient is absent. Figure [5] shows the 
centroid dynamics of an initial Kronccker 5 wave packet 
at a = 3.0 and F = 0.125. The extent of the spatial 
region within which the wave packet asymptotically os- 
cillates is again L' F = W'/F. Thus, we can confidently 
state that Bloch-likc oscillations are rather insensitive to 
initial conditions. 

We studied a biased random tight-binding model where 



the on-site disorder is long-range correlated with the 
power spectrum S(k) ~ l/fc", a > 0. The unbiased 
model supports a phase of delocalized states in the cen- 
ter of the band provided a > 2 Q. We found clear sig- 
natures of Bloch-like oscillations for a > 2 and their ab- 
sence for a < 2. The period of the oscillations agrees 
well with the period in an ideal Bloch band. Thus, the 
amplitude of the oscillations provides a direct way to es- 
timate the energy difference between the two mobility 
edges of the delocalized phase. We found that this en- 
ergy difference is in perfect agreement with the results 
reported in Rcf. This finding opens the possibility to 
experimentally measure the bandwidth of the delocalized 
phase. Actual technological advances made it possible to 
monitor the amplitude of Bloch oscillations in uniform 
semiconductor superlattices Jl9j. Recently, intentionally 
disordered superlattices were used to demonstrate the ab- 
sence of localization in short-range correlated disordered 
systems 0. Thus, we conjecture that intentionally dis- 
ordered superlattices with long-range correlated disorder 
would allow for a clearcut validation of the present re- 
sults. In addition, metallic rings threaded by magnetic 
fields linearly rising in time also display Bloch oscilla- 
tions [l7|. thus opening new experimental ways to test 
our predictions. 
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